NASA Technical Paper 1122 


Simplified Sonic-Boom Prediction 


Harry W. Carlson 

Langley Research Center 
Hampton, Virginia 


ru/\s/\ 

National Aeronautics 
and Space Administration 


Scientific and Technical 
Information Office 




SUMMARY 


A simplified method for the calculation of sonic-boom characteristics for 
a wide variety of supersonic airplane configurations and spacecraft operating 
at altitudes up to 76 km has been developed. Sonic-boom overpressures and sig- 
nature duration may be predicted for the entire affected ground area for vehi- 
cles in level flight or in moderate climbing or descending flight paths. The 
outlined procedure relies to a great extent on the use of charts to provide 
generation and propagation factors for use in relatively simple expressions 
for signature calculation. The computational requirements can be met by hand- 
held scientific calculators, or even by slide rules. With little sacrifice in 
accuracy, complete calculations can often be obtained in less time than is 
required for the preparation of computer input data for the more rigorous cal- 
culation methods. A variety of correlations of predicted and measured sonic- 
boom data for airplanes and spacecraft serve to demonstrate the applicability 
of the simplified method. 


INTRODUCTION 

As the understanding of sonic-boom phenomena has advanced, and the abil- 
ity to provide accurate predictions of sonic-boom phenomena has improved, the 
actual prediction process has become quite complex. For conventional air- 
plane configurations, the usual procedures as described in references 1 and 2 
call for employment of several sophisticated computer programs covering air- 
plane geometry and aerodynamic considerations as well as the wave propagation 
aspects of the problem. For spacecraft operating in the sensible atmosphere, 
it has become the practice (ref. 3) to rely on wind-tunnel tests of small but 
detailed scale models and specialized computer programs for extrapolation of 
the data to full-scale conditions. In either case, the process is complex, 
lengthy, and expensive and requires the services of one or more skilled prac- 
titioners of the art. 

The results of a recent study indicate that for many purposes (including 
the conduct of preliminary engineering studies and the preparation of environ- 
mental impact statements), sonic-boom predictions of sufficient accuracy can 
be obtained by using a simplified method which does not require a wind tunnel 
or elaborate computing equipment. Computational requirements can in fact be 
met by hand-held scientific calculators, or even slide rules. In addition, 
successful use of the method is not highly dependent on the skill and knowl- 
edge of the person performing the calculations. 

This prediction technique results from a simplification of the purely 
theoretical methods described in references 1 and 2, which have been shown to 
provide quite acceptable estimates of sonic-boom phenomena for a wide range of 
flight conditions for conventional airplane configurations. A recent wind- 
tunnel study (ref. 4) has shown that purely theoretical methods may be applied 
to prediction of sonic-boom phenomena for extremely blunt bodies at high 
supersonic speeds, provided that propagation distances are large relative to 



body dimensions. This finding justifies the use of the simplified method for 
spacecraft as well as for conventional airplane configurations. 

This report gives a general description of the prediction method, lists 
the steps required in the calculation procedure, and provides charts that are 
needed in the process. In addition, sample problems are presented to illus- 
trate the method in use, and correlations with flight-test experimental data 
are shown to demonstrate its applicability. 

The reader interested only in prediction of ground-track overpressures 
for conventional airplane configurations may wish to go directly to the sec- 
tion entitled "Further Simplification," in which a brief, self-contained pre- 
diction method is described . 


SYMBOLS 


A(x) area of aircraft cross sections normal to flight direction at a 

given value of x-coordinate (cross sections normal to longitudinal 
axis of aircraft may be substituted in most cases), m 2 


A e (x) 


*e,max 


*e,1 


a v 

B(x) 


®max 

b(x) 


total effective area of aircraft at a given value of x-coordinate, 
A(x) + B(x) , m 2 

maximum effective area, m 2 

total effective area at midpoint of effective aircraft 

p 

length, l ey nr 

speed of sound at aircraft (vehicle) altitude, m/sec 

equivalent cross-sectional area due to lift at a given value of 
x-coordinate, m 2 

p 

maximum equivalent cross-sectional area due to lift, m^ 

local span of aircraft planform at a given value of x-coordinate, m 


d distance between aircraft ground-track position at time of 

sonic-boom generation and location of ground impact point (see 
fig. 5), km 

d x component of d in direction of aircraft ground track, km 

d v component of d in a direction perpendicular to aircraft ground 

track (i.e., in lateral direction), km 


dy c value of dy at lateral limit or cutoff of sonic-boom ground 

* ~ footprint , km 


h altitude of aircraft above ground, h v - h g , km 
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h e effective altitude (see fig. 5), km 

hg altitude of ground above sea level, km 

h v altitude of aircraft (vehicle) above sea level, km 

Kd ray-path distance factor 

K<j,c ray-path distance factor for cutoff conditions, M e = M c 

K<j l00 ray-path distance factor for an infinite Mach number 

Kl lift parameter (see fig. 4) 

Kp pressure amplification factor 

^p,oo pressure amplification factor for an infinite Mach number 

Kp reflection factor, assumed to be 2.0 

Kg aircraft shape factor 

K^ signature duration factor 

K^ joo signature duration factor for an infinite Mach number 

l aircraft characteristic length, normally the fuselage length, m 

l e effective length of aircraft used in determination of aircraft 

shape factor, m 

M aircraft Mach number 

M c aircraft cutoff Mach number below which sonic boom will not reach 

ground 

M e aircraft effective Mach number governing sonic-boom atmosphere 

propagation characteristics 

n<j exponent of Mach number parameter in atmospheric distance factor 

curve fit (see eq. (12)) 

np exponent of Mach number parameter in atmospheric pressure ampli- 

fication factor curve fit (see eq. (13)) 

n^ exponent of Mach number parameter in atmospheric signature duration 

factor curve fit (see eq. (14)) 

p atmospheric pressure, Pa 

Ap incremental pressure due to sonic boom, Pa 
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Apjjgjj incremental pressure at N-wave bow shock, also referred to as bow- 
shock overpressure, Pa 

atmospheric pressure at ground level. Pa 

atmospheric pressure at aircraft (vehicle) altitude, Pa 

aircraft planform area, m 2 

At time increment, sec 

W aircraft weight, kg 

x distance from aircraft nose measured backward along flight path, m 

y flight-path angle (see fig. 5), deg 

0 ray-path azimuth angle (see fig. 5), deg 

<J) angle between aircraft ground track and ground projection of ray 

path (see fig. 5), deg 


METHOD IN GENERAL 

The simplified prediction method is applicable to a wide variety of 
supersonic airplane configurations and spacecraft operating at altitudes up 
to 76 km. It provides estimates of sonic-boom pressure and signature duration 
over the entire exposed ground area for vehicles in level flight or in moder- 
ate climb or descent flight profiles. The effects of flight-path curvature 
and aircraft acceleration, however, are not considered, and the method is 
further restricted to a standard atmosphere without winds. These limitations, 
however, do not appear to affect the general applicability of the method for 
normal variations from the standard atmosphere (temperature profiles and winds) 
and for moderate flight-path curvature and aircraft acceleration. When situa- 
tions do arise in which these effects are important, Whitham F-function data 
provided by the simplified method may be used to supply sonic-boom generation 
data necessary for employment of the propagation computer program described in 
reference 5. 

Another possible restriction to the general applicability of the simpli- 
fied method is the assumption that the pressure signal generated by the air- 
craft is of the far-field type, the classical N-wave (a compression or shock 
wave followed by a linear expansion to pressures below ambient and a second 
shock of equal magnitude which restores ambient pressure) . Normally this is 
a valid assumption, at least for the bow shock and the positive portion of the 
signature, and even when near-field effects become noticeable (some examples 
are shown) the estimate provided by the simplified method tends to be 
conservative. 

The information required for the calculations and the pressure-signature 
predictions provided by the simplified method may be discussed with the aid of 


Pg 

P V 

S 
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figure 1. First, the factors governing the generation of the sonic boom must 
be known. For aircraft not covered by shape-factor charts provided in this 
report, a general description of the aircraft geometry must be provided so 
that shape factors may be calculated. The description need not be detailed, 
and no knowledge of the aircraft lift distribution is required. Of course, 
the aircraft operating conditions - Mach number, altitude, weight, and flight- 
path angle - must be known. In addition, atmospheric pressure and sound speed 
at the aircraft altitude and at ground level must be known. The signature 
data provided by the method include the N-wave bow- shock pressure rise, the 
signature duration, and the location of the ground impact point relative to 
the aircraft position at the time the boom was generated. 

The procedure for calculation of sonic boom by the simplified method 
involves three basic steps: determination of an aircraft shape factor, evalu- 

ation of atmospheric propagation factors, and calculation of signature shock 
strength and duration. These basic steps are implemented as follows: 

( 1 ) Determine aircraft shape factor Kg 

(a) From the aircraft geometry, calculate Kg for a specified 

weight and set of operating conditions according to the 
steps outlined in the section "Shape-Factor Calculation" 

or 

(b) For aircraft covered by the charts discussed in the section 

"Shape-Factor Charts," use aircraft weight and operating 
conditions to calculate lift parameter and read Kg 
directly. 

(2) Determine propagation parameters M e and h e from the operating 
conditions according to equations given in the section "Atmospheric Factor 
Derivation," and read atmospheric factors Kp and Kj. from the charts dis- 
cussed in the section "Atmospheric-Factor Charts." 

(3) Calculate bow-shock overpressure from the equation 

Apmax = K p K R x/^ (m2 - 1 ) 1 /8 h e -3/4i3/4 Ks (1) 


and signature duration from the equation 


3.^2 

At = K t 

a v 


M 

(M2 - 1) 3/8 


he 1/4 z3AK s 


( 2 ) 
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SHAPE-FACTOR CALCULATION 


The computer method of reference 6 has been considerably simplified to 
obtain the following procedure for calculation of aircraft shape factor from 
geometric data. 

The first step, the definition of aircraft cross-sectional areas, requires 
only areas normal to the flight path (except for very large angles of attack, 
areas normal to the aircraft longitudinal axis are equally acceptable) rather 
than those defined by Mach planes. A typical distribution obtained in this 
manner is shown in the following sketch: 



Little loss of accuracy results from this simplification; there is, in fact, 
evidence that for blunt shapes or high supersonic speeds the normal areas are 
more appropriate (ref. 4). Although some care should be taken in establishing 
the maximum cross-sectional area and its location, great precision is not 
required and only the general shape of the remainder of the curve need be of 
concern. If known, the area of the stream tube of air entering the engine 
inlet should be subtracted from the total defined by external contours. 

The second step involves the definition of an equivalent area due to lift. 
Detailed theoretical methods of deriving this contribution are discussed in 
references 1 and 6. It has been found that a reasonably accurate approxima- 
tion of the lift distribution may be defined by the planform area distribution 
as illustrated in the following sketch: 


B(x) 
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The equivalent area due to lift is defined by the equation 


B(x) 



1 W cos y cos 0 
1 . 4p v M 2 S 



b(x) 


dx 


Note that the lifting force which influences the sonic boom is defined by an 
aircraft weight component normal to the flight path and directed along the 
initial ray-path azimuth angle 0 . 

The third step is the combination of these two area contributions to 
obtain a distribution for the total effective area of the aircraft, illus- 
trated in the following sketch: 



From this curve, the maximum effective area A e max and its location can be 
determined. Then, as shown in the sketch, the effective length l e and the 
effective area A e>1 required for evaluation of the aircraft shape factor may 
be selected. 

Finally, in the fourth step, the aircraft shape factor may be found from 
a reading of the shape-factor parameter curve of figure 2 and insertion of the 
appropriate areas and lengths. The shape-factor parameter has been defined by 
evaluation of sonic-boom theory for an effective area distribution of the form 
A e (x) = k-jx + k£X 2 , with constants k-] and k£ selected so the curve passes 
through the A e 1 and A e(max points. (See sketch at top of fig. 2.) Note 
that the aircraft shape factor is primarily dependent on the maximum effective 
area. The location of the maximum area and the shape of the effective area 
development play lesser roles. 

The ability of the simplified procedure to account for shape-factor vari- 
ation of a representative supersonic airplane is illustrated in figure 3. 
Effective area developments shown in the figure cover the extremes and the 
midpoint of lift parameters likely to be encountered in normal flight. Signa- 
tures shown below each of the effective area developments indicate the influ- 
ence of the shape-factor variation on the signature characteristics. As can 
be seen in equations (1) and (2), both overpressure and duration vary directly 
with the magnitude of the shape factor. It has been found that the simplified 
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curve-fit method gives shape factors which are generally within 5 percent (and 
almost always within 10 percent) of the values given by the more rigorous theo 
retically based computer methods. 


SHAPE-FACTOR CHARTS 

For the reader's convenience, aircraft shape factors for a variety of air- 
planes employed in sonic-boom flight-test programs and for more contemporary 
aircraft are given in the charts of figure 4. To use these charts it is nec- 
essary only to identify the aircraft of concern, calculate the lift parameter 
K L , and read' the shape factor directly. Note that with the exception of the 
space shuttle orbiter, which is set apart as a result of its high volume, shape 
factors for all the aircraft fall in a fairly narrow band. For aircraft not 
specifically covered by the charts, shape factors may be chosen by selecting 
a similar configuration. Generally, larger airplanes are more slender and will 
have lower shape factors. Fighter airplanes, on the other hand, tend to have 
higher shape factors. At the upper end of the lift-parameter range, the shape 
factor is more dependent on the wing planform than on the aircraft size and 
volume, and highly swept wings tend to have better shape factors. 

There is no great sensitivity of the final overpressure and duration cal- 
culations to the choice of the aircraft characteristic length l, even though 
it appears as a squared term in the lift parameter. It is important, however, 
that the chosen length be used consistently throughout the calculations. 

If, as mentioned in the introduction, it becomes necessary to consider 
the effects of curved flight paths and acceleration in situations where far- 
field conditions are to be expected, the shape factor provided by this simpli- 
fied prediction method may be used to supply an input F-function for the atmo- 
spheric propagation program of reference 5. The necessary F-function may be 
represented by a double triangular distribution as shown in the following 
sketch: 


F-function 
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Definition of the lifting force normal to the flight path used in shape-factor 
determination must now include the effects of aircraft accelerations; it is no 
longer appropriate to use only a component of the aircraft weight as is done 
throughout this report. 


ATMOSPHERIC-FACTOR DERIVATION 

Atmospheric propagation factors for ray-path distance, pressure amplifica- 
tion, and signature duration have been derived from repeated use of the sonic- 
boom propagation computer program described in reference 5 for a matrix of 
altitudes, Mach numbers, ray-path azimuth angles, and flight-path angles. To 
determine these factors, the program was run first with a standard atmosphere, 
and then with a uniform atmosphere with a pressure equal to the geometric mean 
of the s tanda rd atmospheric pressures of the aircraft altitude and at ground 
level, \J p v Pg. The factors are simply the ratios of sonic-boom parameters 
given by the first run to those given by the second. 

In order to broaden the applicability of the atmospheric propagation 
factors and allow the use of on-track and level-flight data for off-track and 
flight-path-angle conditions, the concepts of effective Mach number and effec- 
tive altitude were employed. It was presumed that the ray-path curvature and 
the overpressure amplification due to atmospheric effects would depend primar- 
ily on the initial inclination of the aircraft-generated ray path with respect 
to the earth plane and on the ray -path distance measured perpendicular to the 
aircraft flight path. A convenient measure of the ray-path inclination is the 
Mach number for level flight, which would have the same ray-path angle in the 
flight-track plane. This is termed the effective Mach number. The ray-path 
distance measured perpendicular to the aircraft flight path is termed the 
effective altitude. Geometric relationships used in the derivation of these 
parameters are illustrated in figure 5. Simplified equations for effective 
Mach number, propagation distance, and effective altitude for the on-track 
case are as follows : 


M e = 


sin 


(y + cot -1 v/mTH) 


(3) 


<3 X = K d| 


vK 2 - 


(K d from fig. 7(b)) 


(4) 


h e = h cos y + d x sin y 


(5) 
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Complete equations for the more general off-track case are as follows: 


d = 



<{) = tan 


1 tan 0 cos y ( 1 + tan2 Y) 


tan y + 


cos 0 y/52 - 1 


d x = d cos <t> 


dy = d sin 


h e = sj d y 2 + (h cos y + d x sin y)‘ 


(8) 

(9) 

( 10 ) 

(ID 


The usefulness of the effective Mach number and effective altitude in 
reducing a large amount of computer-generated data to manageable proportions 
xs illustrated in figure 6. For a selected altitude of 20.0 km, ray— path dis- 
tance factor and pressure amplification factor Kp are shown a3 a func- 

tion of the effective Mach number. The Mach number scale was chosen in order 
to cover Mach numbers from one to infinity with an expansion of the lower Mach 
number range. Program-generated data are given for ray-path azimuth angles 
from 0° to 60° and for flight-path angles from -15° to 30°. The nesting of 
the data verifies the concepts of the effective Mach number and the effective 
altitude. For other altitudes up to about 70 km, the tendency toward collapse 
of the data to a single line is generally somewhat more pronounced than that 
shown in figure 6. 


ATMOSPHERIC-FACTOR CHARTS 

Atmospheric factors required in this simplified method for sonic— boom 
prediction are presented in figure 7. The first chart (fig. 7(a)) gives the 
cutoff Mach number in a standard atmosphere as a function of aircraft altitude. 
Cutoff occurs for disturbances which propagate away from the aircraft along a 
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ray path (as defined by the effective Mach number) which curves under the 
influence of atmospheric gradients to an extent just sufficient for it to 
become horizontal at ground level. As is discussed later, this cutoff limits 
the lateral extent of the ground area affected by the sonic boom. If the 
effective Mach number calculated according to the appropriate equation from 
the preceding section is not greater than the cutoff Mach number for the air- 
craft altitude, there is no need to go further; the signal will not reach the 
ground . 

Ray-path distance factors are given in figure 7(b). The horizontal 
distance traveled by a particular ray path from the time a signal leaves the 
aircraft until it reaches the ground may be found by using this factor. The 
factor read from the chart for a given altitude and effective Mach number 
is simply multiplied by the distance traveled in a uniform atmosphere 

(h /\j Mg^ - 1 j . Interpolation must be done with care because of the nonlinear 
nature of the data. There is no need to use the ray -path distance-factor 
charts for estimates of level- flight on-track overpressure; they are required, 
however, to obtain an effective altitude for overpressure calculation at off- 
track locations and for flight-path angles other than zero. 

The limiting or cutoff ray-path distance factor may be read from fig- 
ure 7(c). The primary use of this chart is in defining the lateral extent 
dy q of the affected ground area for level flight given by 


d y.c 




2 


The chart is also useful in defining the sonic-boom footprint limits for mod- 
erate climbing or descending flight profiles. These limits, however, must be 
found by iteration. The cutoff distance chart was prepared from an analysis 
given in reference 7- The computer program of reference 5 does not provide 
this information since calculations are terminated when ray paths come within 
2° of the horizontal. 

Pressure amplification factors are presented in figure 7(d). Care is 
required for accurate interpolation of the charts because of the nonlinear 
nature of the data. This factor, in combination with the effective-altitude 
term, accounts for nonuniform atmospheric effects in the equation for bow- 
shock overpressure given in equation (1). 

The final atmospheric-factor chart (fig. 7(e)) gives the signature dura- 
tion factor which, in combination with the effective-altitude term, accounts 
for the effects of a nonuniform atmosphere in the equation for signature dura- 
tion given in equation (2). The actual Mach number rather than the effective 
Mach number is used when reading this chart. 

The atmospheric-factor charts presented herein apply directly to the sit- 
uation where ground level is at or near sea level. The charts can, however, 
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be used with little error for ground levels up to about 1600 meters. The air- 
craft altitude above sea level h v should be used when entering the charts, 
and the altitude above ground level h should be used when calculating the 
effective altitude and propagation distances. 

Generally, the atmospheric propagation factors may be read directly from 
the previously discussed charts. However, the charts which cover only Mach 
numbers greater than 1.2 do not give sufficient information for situations in 
which cutoff conditions are approached. To cover these situations, the follow- 
ing curve-fit equations which employ exponents given in figure 8 may be used: 


*d 


Kd,c + (Kd,„ 




( 12 ) 



(13) 


Kt 



(14) 


The form of these curve-fit equations was suggested by the nature of the pro- 
gram data as exemplified in figure 6. The asymptotic approach of the factors 
to limiting values at Mach numbers approaching infinity should be noted. The 
pressure factor displays singular behavior as the cutoff Mach number is 
approached, but the distance factor has a finite limit at this Mach number. 
These characteristics were taken into account in the formulation of the curve 
fit. An example of the curve-fit approximation is given in figure 6. 

Equations (12) to (14) and the exponent curves of figure 8 appear to be 
somewhat formidable. The information, however, may be employed in a rather 
straightforward manner. First, it is necessary to define the cutoff Mach num- 
ber and limiting values of the atmospheric factors from figure 7. Then, for 
the given altitude, read the curve-fit exponents for each of the three factors 
from the chart of figure 8. Finally, find the atmospheric factors by substi- 
tuting the desired effective Mach number, the cutoff Mach number, the limiting 
factors, and the curve-fit exponents into equations (12) to (14). 


SIGNATURE CALCULATION 

In the preceding sections, means of obtaining an aircraft shape factor 
were described and the procedure for reading atmospheric propagation factors 
from the charts was outlined. In the process of evaluating the atmospheric 
propagation factor, the location of the impact point on the ground relative to 
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the position of the aircraft at the time of generation (d x and dy) is also 
provided. 

With the shape factors and atmospheric factors in hand, the signature may 
be calculated by applying equations (1) and (2). All the terms are described 
in the section "Symbols" and some are discussed in detail elsewhere in this 
report . 

It may be of help here to discuss the reflection factor Kp. When the 
shock wave meets the Earth, flow velocities behind the shock are altered; 
vertical velocities are blocked and their energy is converted into a pressure 
increase. For weak shocks meeting a smooth, rigid, and flat surface, the 
resultant overpressure is theoretically twice that of the free-air value of 
the incident wave. In flight tests, reflection factors varying from about 1.8 
to 2.0 have been recorded. The lower values are generally associated with 
sandy or marshy terrain. For ground overpressure predictions presented in 
this report, a reflection factor of 2.0 has been assumed. 

Although the calculative procedures are relatively straightforward and 
simple, it may be desirable to utilize the programming capabilities of some of 
the modern pocket-sized or desk-top calculators for repetitive calculations. 
Some notes on programming are offered in appendix A. 

A series of three sample problems has been prepared to aid the reader in 
understanding the calculation process. These problems are presented in 
appendix B. 


CORRELATION WITH FLIGHT -TEST DATA 

To demonstrate the applicability of the simplified sonic-boom prediction 
method to a variety of problems, correlations with flight-test data (refs. 8 
to I 1 !) are shown in figures 9 to 13. 

In figure 9, predictions are compared with measured signatures for three 
distinct airplane types (refs. 8 and 9) operating at or near cruise conditions. 
Only for the fighter aircraft, which displays some evidence of near-field 
effects, is there an appreciable difference between predicted and experi- 
mental results. In this situation, as noted previously, the prediction is 
conservative. 

Correlations of predictions with measurements for a reconnaissance air- 
plane (ref. 10) covering a range of flight altitudes are shown in figure 10. 

The variation of Ap^y with h v is quite typical of that for most airplanes; 
the Mach number differences are of secondary importance. 

The validity of the method for off-track predictions for a bomber air- 
plane (ref. 8) is illustrated in figure 11. At large lateral distances, the 
ray paths curve sufficiently to approach the Earth at near-grazing angles. 

This results in signature distortions and ultimately in a breakup of a recog- 
nizable signature into random noise beyond the lateral cutoff. There is some 
evidence of such a transition for the most distant signature shown in the 
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figure. In general, the variations with lateral distance are predicted reason- 
ably well; however, there is a singularity predicted at the cutoff point which 
does not appear to materialize. Other examples of measured lateral-spread 
patterns are given in reference 15. No evidence of pressure amplification at 
the cutoff point is shown in reference 15, and the authors apparently have 
chosen to fair the prediction curve in a monotonically decreasing fashion 
rather than attempt to define the theoretical singularity. Users of this 
system may choose to do likewise. 

It is remarkable that the predictions for the Apollo command modules 
shown in figure 12 bear any similarity to the measured data (refs. 11 and 12). 
The diameter-length ratio of the effective body shape is about 9 and the 
effective length is only 0.42 m. The appearance of the unexpected multiple 
shocks in the measured signatures remains unexplained. Although the correla- 
tions are far from exact, the predictions are nevertheless useful. Especially 
notable is the degree of correlation indicated for a Mach number of 15. 

Another example of the applicability of the simplified prediction method 
to spacecraft is given in figure 13. As mentioned in the discussion of the 
sample problems (appendix B) , the sonic boom generated by the Apollo launch 
vehicles (refs. 12 and 13) is dominated by the exhaust gas plume, so this is 
a very special situation. Another special feature is the ground pattern of 
the signatures. As shown in the flight-path profile there are two families 
of rays. One family consists of near-grazing rays generated prior to genera- 
tion of that ray which reaches a point nearest the launch. The other family 
consists of more direct rays generated after that time. As can be seen on 
the plot of overpressure versus distance from launch, there is a portion of 
the ground track over which double or superimposed signatures are predicted. 

It is apparent from the inset pressure-signature sketches that the grazing- 
ray component of these signatures is grossly overpredicted. However, the 
pressure-doubling at the closest ground intersection point appears to have 
materialized. Large distances traveled in proximity to the Earth appear to 
cause disintegration of a distinct signature into random noise. There is 
also a question as to whether grazing rays should be considered to have a sur- 
face reflection factor of 2 or a free-air reflection factor of 1 . These data 
seem to indicate that the grazing-ray pressure signatures should be ignored. 
However, at the first ground-track point, members of the grazing-ray family 
have not yet reached grazing conditions and thus contribute to the signature. 
With this stipulation, the predicted overpressure indicated by the solid line 
agrees reasonably well with the measured data. The predicted signature dura- 
tion, however, as indicated by the sketches, is too large by a factor of 
nearly 2. The curved flight path and the vehicle acceleration may play some 
part in the cause of the discrepancy. 


FURTHER SIMPLIFICATION 

For readers interested only in prediction of on-track overpressure for 
conventional airplane configurations, several of the steps in the general 
simplified method may be omitted. In this case it is possible to place all 
the necessary information, with the exception of atmospheric tables, on a 
single page. This "super simplified" method is given in figure 14. 
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Predictions obtained by this method are compared with flight-test data for 
three airplanes (refs. 8 and 14) in figure 15. Each prediction is shown as a 
band to account for variations in airplane weight and Mach number that occur 
at a given altitude. Deviations of predictions from measurements occur at the 
lower altitudes for both the large bomber and the small fighter. These dif- 
ferences arise because the simplified methods do not account for the presence 
of a small amount of near-field effect in the signatures generated by these 
airplanes. Such discrepancies are anticipated for very large and slender air- 
planes, but the measured signatures for the small fighter are not at all typi- 
cal of fighter aircraft. The near-field effect for the fighter is due to the 
nontypical long slender nose. A redeeming quality of the prediction method 
is that when near-field effects do occur, the prediction tends to be 
conservative . 


CONCLUDING REMARKS 

A simplified method for calculation of the sonic-boom characteristics of 
greatest interest for a wide variety of airplane configurations and spacecraft 
has been described. The procedure, which has been outlined in a step-by-step 
manner, relies to a great extent on the use of charts to provide the necessary 
sonic-boom generation and propagation factors for use in relatively simple 
expressions for signature characteristics. Computational requirements can be 
met by hand-held scientific calculators or by slide rules. With little sacri- 
fice in accuracy, complete calculations can often be obtained in less time 
than is required for the preparation of computer input data for the more 
rigorous calculation methods. A variety of correlations of predicted and 
measured sonic-boom data for airplanes and spacecraft serve to demonstrate 
the applicability of the method. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 31, 1978 
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APPENDIX A 


NOTES ON PROGRAMMING 

As was mentioned in the main body of the report, it was found that the 
computational needs of the simplified sonic-boom prediction method could be 
met by hand-held scientific calculators. With programmable pocket-sized or 
desk-top calculators, the labor of performing the calculations can be reduced 
to an almost negligible level. Some notes on strategy for employment of pro- 
grammable calculators are offered here. 

Before a programmed solution may be employed, it is necessary to have 
shape-factor charts for the aircraft of concern. For aircraft not covered in 
figure 4 , a shape-factor chart may be prepared according to the steps outlined 
in the section "Shape-Factor Calculation." 

It was found to be convenient to calculate the atmospheric factors first. 
The primary input quantities are the flight conditions M, 9, Y> an d h v 
and the ground level h g . If the curve-fit data are to be used in determina- 
tion of the atmospheric factors, the quantities Mq, Kd,oo' ^d.q* ^p,a>» 

K* ru nr,, and n* read from the charts for the flight altitude may also 

be'Siedts initial inpSt. 

The effective Mach number M e and the ray -path azimuth angle <t> may be 
calculated as follows: 


M e = 


1 + 


[A(1 - B tan y)j 2 
[A(tan y + B )J 2 + (CD) 2 


<{> = tan - ^ 


^tan 9 cos y D' 
^ tan y + k , 


where 


1 


A = 


\|m 2 - 


cos y \JM<- - 1 
1 


B = 


C = 


cos 9 \j M 2 — 1 
tan 9 

n|m 2 - i 


D = tan 2 y + 1 
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With the effective Mach number known, the atmospheric factors may be determined 
either from a direct reading of the charts (which requires a program temporary 
stop for chart reading and input of K d , K p , and K t ) or from the curve-fit 
data. When evaluating atmospheric factors from the curve-fit data, first test 
M e to see if it is greater than M c . If it is not, the calculation must be 
terminated; the ray path chosen will not reach the ground. To evaluate the 
atmospheric factors K d , FC p , and from the curve-fit data, use the follow- 

ing expressions: 


/M _ M \ n d 

^d = K d , c + (K d(m -K d>c )^-f_ij 

/M e - 1 \ n P 
Kp = ) 

Kt = Mm— ) 


The ground location of the boom and the effective altitude may now be 
found as follows: 


K d (h v - h K ) 

d = 

\K 2 - ' 

d x = d cos <{> 


dy = d sin <J> 


h 


e 



hg) cos y + d x sin 


To calculate the magnitude of the sonic-boom overpressure, first evaluate 
the lift parameter as follows: 


KL 



- 1 W cos y cos 6 
1.Hp v M 2 ! 2 


When Kl is known , it is necessary to provide another temporary stop for 
reading and input of the aircraft shape factor from the appropriate curve. In 
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addition, the quantities p v , Pg, a v , and Kr must be input. Now the bow- 
shock overpressure and the signature duration may be calculated as follows : 

APmax = K p K R \/PvPg < M 2 “ 1 ) 1 /8 h e " 3/1 h 3/4 Ks 


At = K t h e 1/4z3/4 Ks 

a v ( M 2 _ i) 3/8 


Any consistent set of measurement units may be used as input. The over- 
pressure will be in the units of the atmospheric pressure and the signature 
duration will be in the time unit used for the speed of sound. 
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SAMPLE PROBLEMS 

The calculations required for the three sample problems are presented in 
this appendix with only a brief explanation. The first sample problem illus- 
trates the calculation of the shape factor from aircraft geometric data, in 
this case taken from small-scale three-view drawings such as those found in 
reference 16. The calculations shown in the second problem cover the situa- 
tion in which an aircraft shape factor may be read from the charts. The off- 
track signature calculation, however, is somewhat more complicated than the 
on-track calculation of the first sample problem. The third problem covers 
an unusual situation in which the sonic-boom behavior is dominated by the 
exhaust gas plume of the propulsion system rather than by the vehicle itself. 
This problem also demonstrates the versatility of the simplified method. The 
area of the exhaust gas plume of the ascending Apollo launch vehicle was 
obtained from data presented in reference 17. 


Sample Problem 1 : Reconnaissance Airplanes, On-Track 

The flight conditions and geometric data are 

M = 1.99 W = 41 000 kgf = 402 000 N 

1 = 32.7 m Y = 0° 0=0° 


h v = 15.240 km p v = 11.6 kPa a v = 295 m/sec 

h g = 0.760 km p g = 92.3 kPa K R = 2.0 

h = h v - hg = 14.480 km 

The aircraft shape factor is obtained as follows: 


Jm 2 - 1 W _ 

Bmax = ~ ~ 10.7 m 

1 .4p v M 2 
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From the area plots at right 
(from data given in ref. 16) 


APPENDIX B 


l e = 24 m 
A e,max = 15.0 m 2 

A e ,1 = 4 - 4 m 2 

A e , 1 / A e , max = 0 • 2 ^ 


From figure 2, 

K S = 0.685 


\J A e , max 
z3/4 z 1/4 



= 0.087 


Propagation factors are obtained from equations (3) to (5) and figure 7. 
Me a M * 1.99 
h e = h = 14.480 km 
K p = 1.10 
K t = 0.85 

The signature is calculated as follows: 

APmax = K p K R \/PvPg (m2 - 1 ) 1 /8 h e “ 3/4 i 3/4 Ks = 74 Pa 

At = Kt — h e 1/4z3/4 Ks = 0.17 sec 

a v (m 2 - 1) 3/8 
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Sample Problem 2: Medium Bomber Airplane, Off-Track 

The flight conditions and geometric data are 

M = 2.0 W = 36 000 kgf = 353 000 N 

l = 30 m Y = 0° 0 = 52° 

h v = 18.600 km p v = 6.88 kPa a v = 296 m/sec 


h g = 0.760 km p g = 92.3 kPa K R = 2.0 

h = h v - h g = 17.840 km 

The aircraft shape factor is obtained from figure 4. 


- 1 W cos y cos 0 
1 .4p v M 2 ! 2 


Kg = 0.088 


Propagation factors are obtained from equations (6) to (14) and 
figures 7 and 8. 



dy = d sin <)> = 30.2 km 
h e = Jd y 2 + h 2 = 35.1 km 
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M c = 1.153 

/Me - M o\" d 

Kd=(^— rj < K d,„- K d , 0 ) - K<J,0 = '.32 



1 .42 


Kt = K t 



0.82 


The signature is calculated as follows: 


APraax = Kp K R\fp^ - 1 ) 1 /8 h e - 3/4 i 3/4 K s = 36.1 Pa 

At = K t h e 1/2l i3/4 Ks _ o.19 sec 

a v (M 2 - 1) 3/8 
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Sample Problem 3: Apollo Launch Vehicle in Ascent 

The flight conditions and geometric data are 


M = 4.08 
l = 104.5 m 
h v = h = 36.1 km 
h g = 0 


125 sec and 32.9 
Y = 31 .6° 
p v = 0.490 kPa 
Pg = 101 kPa 


km from launch 
9=0° 

a v = 310 m/sec 

K r = 2.0 


The aircraft shape factor is deter- 
mined from area plot at right (from data 
given in ref. 17). 



From figure 2, 


K S 


0.65 


yp 


‘e,max 


z3/4 Ze 1/4 


0.68 


/ 

,max 
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Propagation factors are obtained from equations (3) to (5) and figure 7. 

M e = 7 1 . — v = 1.39 K d = 0.99 

sin (y + cot _ 1 v - ij 

d x = K d f , h ] = 37.2 km K p = 1.40 

h e = h cos y + 9 X sin y = 50.2 km = 0.73 

Impact point = 32.9 + d x = 69.3 km from launch 
The signature is calculated as follows: 

A Pmax = K p K R \fPvPg (M 2 * 1 ) 1 /8 h e - 3/4 i 3/4 K S = 184 N/m 2 


At = K t 


3.42 M 

3v (M 2 - 1) 3/8 


h e 1/ V / ' , Ks 


3-90 sec 


The reader may be concerned that the effective-area curve terminates at 
a point where the area is still growing. This, however, does not cause any 
real difficulty. The major compression of the airflow about this shape takes 
place in the first two-thirds of the chosen effective body length. Beyond 
that point only relatively mild compressions or expansions occur. For example, 
if the shape given in reference 9 were extended along the dashed line (which 
appears to be a reasonable assumption) to an effective length 50 percent 
greater than the original effective length, the signature overpressure and 
duration would be unchanged. Also, the chosen effective length l e does not 
in this case begin at the origin of the area plot. 
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Basic information required 


Aircraft geometry (for aircraft not covered 
by shape-factor charts) 

Operating conditions 
Mach number 
Altitude 
Weight 

Flight -path angle 

Atmospheric data 

Ambient pressure 
Sound speed 



Figure 1.- Information required and data provided by simplified 
sonic-boom prediction method. 
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Curve fit 



Figure 2.- Shape-factor parameter. 
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Time increment, At Time increment, At Time increment, At 

Figure 3.- Comparison of simplified procedure approximation with computer program 

results for actual effective areas. 
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Figure 4.- Representative aircraft shape factors. 




0.74 JKr + 0.027 
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Figure 4.- Concluded 
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(b) Ray-path distance factor 
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(c) Cutoff ray-path distance factor 







Lification factor 
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(e) Signature duration factor 
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Figure 8.- Atmospheric-factor curve-fit exponents. 










Experiment Prediction 



Figure 9-- Correlation of predictions with measurements covering a range 

of airplane types. 
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Figure 10.- Correlation of predictions with measurements covering a range 
of altitude for a reconnaissance airplane. 




Experiment 



H3 


Figure 11— Correlation of predictions with measurements covering a range 
of lateral distances for a bomber airplane. M = 2.0, h„ = 18.5 km 





Figure 12.- Correlation of predictions with measurements made during 

reentries of Apollo 15 and 16. 
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Figure 13.- Correlation of predictions with measurements made during 

ascents of Apollo 16 and 17. 
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Figure 14.- Super-simplified sonic-boom prediction method for on-track 
bow-shock overpressure of conventional airplanes in level flight. 
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Figure 15.- Correlation of on-track bow-shock overpressure predictions 

with measurements. 
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